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Abstract 


We study theoretically potential advantages of two-component mixtures in optical lattices with state-dependent tunneling for ap¬ 
proaching long-range-order phases and detecting easy-axis antiferromagnetic correlations. While we do not find additional advan¬ 
tages of mixtures with large hopping imbalance for approaching quantum magnetism in a harmonic trap, it is shown that a nonzero 
difference in hopping amplitudes remains highly important for a proper symmetry breaking in the pseudospin space for the single- 
site-resolution imaging and can be advantageously used for a significant increase of the signal-to-noise ratio and thus detecting 
long-range easy-axis antiferromagnetic correlations in the corresponding experiments. 
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1. Introduction 

Due to a recent experimental realization of state-dependent 
optical lattices for two-component mixtures of ultracold 
atoms with the magnetic-field-gradient technique and low heat¬ 
ing rate [1] it is now much easier to access and study asymmet¬ 
ric lattice models without a necessity of using heteronuclear 
fermionic mixtures (e.g., ®Li-"^°K) or long-living metastable 
electronic states of the same fermionic isotope (e.g., 37*0 state 
of '^^Yb). Among potential applications of this technique 
one can suggest approaching long-range magnetically-ordered 
states [2]. It is known that quantum magnetism in ultracold 
fermionic mixtures is one of major experimental challenges 
nowadays and a significant progress already has been made 
in this direction. In particular, short-range antiferromagnetic 
(AFM) correlations were effectively measured [3, 4] and their 
unique dynamics in the presence of the tunable lattice geometry 
was observed recently [5]. 

Considering two-component fermionic mixtures from the 
point of view of theoretical models and spin symmetries, opti¬ 
cal lattices with state-dependent (i.e., spin-dependent) hopping 
amplitudes effectively break the initial continuous SU(2) sym¬ 
metry of the system described by the Hubbard model towards 
U(l)xZ 2 , where I 2 is a discrete reflection symmetry along the 
easy axis. The easy-axis direction is important, in particular, 
for the experimental detection of AFM correlations based on 
the Bragg spectroscopy analysis [4, 6] and on the quantum-gas- 
microscope technique (QGMT). Despite the fact that tempera¬ 
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tures and entropies* achieved with recent successful develop¬ 
ments of the QGMT for fermionic mixtures [7, 8, 9, 10, 11, 12] 
are high to observe long-range magnetic correlations, thus fur¬ 
ther optimizations and improvements in cooling protocols are 
required, it is important to study characteristic dependencies 
of these thermodynamic quantities on other system parame¬ 
ters (including different symmetries of magnetic ground states), 
thus determine the most optimal regime for in situ imaging of 
the long-range AFM correlations. 

2. Theoretical description 

Ultracold two-component fermionic mixtures in optical lat¬ 
tices with a sufficiently strong lattice potential, V\^t > SE,, 
where E, is the recoil energy of atoms, are well described by 
the single-band Hubbard model with the Hamiltonian 

77= “ 2j + h.c.) + U X «iT”a 

<ij) o- i 

+yYj^nlafni-Y^Y^fio-nia; ( 1 ) 

i i (T 

where to- is the hopping amplitude of fermionic species in a 
particular atomic hyperfine state that we denote by the spin-1/2 
index cr = {T^sl), thus we consider here and below two atomic 
components as pseudospins, (c/o-) is the corresponding cre¬ 
ation (annihilation) operator of atoms at the lattice site /, the 


^ The entropy per particle is more crucial quantity in the context of ultracold- 
atom experiments, since the system does not exchange heat with environment 
and its parameters can be changed adiabatically. 
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notation {ij) indicates a summation over nearest-neighbor sites, 
and U is the magnitude of the on-site repulsive (U > 0) inter¬ 
action of the two different species with corresponding densities 
n,T and hji {hia- - c,o-). In the third term, V is the amplitude 

of the external harmonic potential, r, is the distance from the 
lattice site i to the trap center, and a is the lattice spacing. In the 
last term, is the chemical potential that determines the total 
number of atoms of each spin component in the system. 

Note that the asymmetric Hubbard model (1) (it can also be 
recognized as the extended Falicov-Kimball model for spin¬ 
less fermions in the context of solid-state materials [13]) can be 
transformed to an anisotropic Heisenberg (or a spin-1/2 XXZ) 
model in the limit of Ujta- » 1 and tii ^ 1. The latter is de¬ 
scribed by the Hamiltonian 

(ij) iU) 

( 2 ) 

i 

with the constants /y = 2(f^ + Ji. - Ap = 

(//| — and the spin-1/2 operators Sf - (here and 

below we use units h - 1), where cr* are the Pauli matrices 
{R = {X,Y,Z}). Hence, we see that the presence of hopping 
imbalance (f| results in /y > 7^, thus breaks the SU(2) 
rotational spin symmetry towards U(l)xZ 2 , where the discrete 
symmetry Z 2 can be broken either spontaneously by long-range 
AFM ordering along the Z axis or by the chemical potential dif¬ 
ference A// 0 that plays a role of the external magnetic field 

favoring the ferromagnetic configuration along the same axis. 
Naturally, with an increase of A/j. at the fixed asymmetry in hop¬ 
ping amplitudes the easy-axis (Ising-type) AFM configuration 
becomes less and less energetically favorable, thus a transition 
to another AFM-ordered easy-plane (canted) many-body state 
becomes possible (see Refs. [14, 15, 16] for more details). 

Therefore, to avoid a potential competition between differ¬ 
ent types of AFM ordering that can also result in a significant 
suppression of critical temperatures, in the original model (1) 
we consider chemical potentials the same for both spin com¬ 
ponents (i.e., jji - fii = fi). Note that asymmetric hopping 
amplitudes fy ti together with jUy = result in a nonzero po¬ 
larization, i.e. not equal total numbers of particles in two spin 
states (Ay + Ay), of the trapped system [2, 16], but this condi¬ 
tion is the most optimal for the easy-axis AFM ground state of 
the model (1) at half filling (e.g., at yU = t//2 and r,- = 0), as 
discussed above. 

Below, we consider a three-dimensional optical lattice setup 
with the Hubbard parameters entering Eq. (1) that are set close 
to the experimental values [1]. In particular, we focus on two 
opposite limits: (i) zero (or very small) hopping imbalance fy = 
fy = f and (ii) large hopping imbalance (e.g., fy = 0.54f and 
fy = 0.06f)- According to Ref. [1], both cases can be effectively 
realized by the magnetic-field-gradient technique with a high 
level of control. 

Our theoretical analysis is based on the dynamical mean-field 
theory (DMFT) [17] with the exact diagonalization solver [18] 
and the number of orbitals Hs - 5 per each spin component 


in the corresponding Anderson impurity model. Since we are 
interested mostly in the easy-axis observables, in DMFT it is 
enough to account for the standard hybridization terms between 
the impurity and the bath [18]. The corresponding Anderson 
parameters of the impurity model are found iteratively till the 
convergence based on DMFT self-consistency conditions [17] 
is reached. For the given values of the Hubbard parameters 
our approach allows to calculate the local observables, such as 
expectation values of the density of particles of any spin com¬ 
ponent, the double and the hole occupancy, as well as the fluc¬ 
tuations of the particle number on a particular lattice site. 

To account for the inhomogeneity effects produced by the 
external trap, we use DMFT with the local density approxima¬ 
tion (LDAh-DMFT). Note that EDA does not account for the 
proximity effects close to the phase boundaries, however in the 
cases under study these effects do not play a crucial role lead¬ 
ing only to minor corrections. Within EDA we obtain the local 
observables at the specific point r of the trap from the condition 
- Ilf) - V{rfa)^, where fio is the chemical potential in the 
trap center that for the fixed values of the Hubbard parameters 
defines also the total number of particles in the system. 

From the converged solutions of LDA-hDMFT on different 
lattice sites (i.e., with different r) one can analyze the de¬ 
pendence of the local observables on the distance r (see also 
Ref. [19] for more details). In particular, by combining the re¬ 
sults with the Maxwell relation dsld/u - dnIdT we obtain the 
entropy per lattice site at the particular point ro of the trap (for 
simplicity, we use the units of = 1 and a - I below) 


i(ro) = 2y 



dn(r, U, T) 


rdr. 


(3) 


where the cut-off distance Rmax is determined from the condi¬ 
tion U, T) - 0. 

A subsequent integration of the entropy and density distri¬ 
butions in the trap determine the total entropy S and the total 
number of particles A in the system (here and below we assume 
the axial-symmetric three-dimensional setup) 

^^max 

s(r)47rr^dr, A = I nirfAjiP'dr. (4) 
0 Jo 


Both quantities, S and A, can be considered as the preserved 
numbers in the experiment (and, in particular, during the lattice 
ramp) that allows to access the initial values for the entropy and 
temperature. Note that, alternatively, one can also introduce an 
additional term corresponding to the amount of entropy per par¬ 
ticle As that is added to the system due to uncontrolled heating 
processes during the lattice ramp, as it was done in Ref. [3]. 
Below, for simplicity and consistency reasons, we consider that 
the change in system parameters can be performed adiabatically 
(As = 0). 

In order to determine the initial temperature T in the system 
that is necessary for observations of the many-body quantum 
phases under study we use the expression for the entropy of the 
Fermi gas under assumption of a moderate scattering length a, 
(kplfl,! < 1/2) [20] 

S ^ Ntt^TITp. (5) 
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Figure 1: Dependence of the entropy per particle (left axis) and the initial 
temperature (under assumption of As — » 0; right axis) on the temperature in 
a simple cubic lattice obtained by DMFT for two-component mixtures with 
two different realizations of imbalance in hopping amplitudes, U = lOf*, and 
V = 0.015t. The gray-shaded bar in the bottom paif qualitatively indicates the 
strength of the AFM correlations in the regions with « » 1. The reference value 
of for the balanced mixture is taken from Ref [21], 

Eqs. (3)-(5) allow to set a direct correspondence between ther¬ 
modynamic quantities before and after the lattice ramp. There¬ 
fore, the problem can be effectively solved under assumption of 
adiabaticity of the ramp process. 

3. Results 

3.1. Comparative entropy analysis of the trapped mixtures in¬ 
cluding hopping imbalance 

To proceed with the quantitative analysis, let us specify the 
Hubbard parameters, which we set maximally close to the ex¬ 
perimental values [ 1 ]. In particular, to simplify the theoretical 
description and in accordance with Eqs. (3)-(5), we consider 
the axial-symmetric three-dimensional trap with the frequency 
d) - 2n X 68.4 Hz (i.e., the geometric mean of three frequen¬ 
cies cui, 2,3 in [!])■ By specifying t = 174 Hz [1], we obtain 
the amplitude of the trapping potential V - 0.015f in Eq. (1). 
Concerning imbalances in hopping amplitudes, we focus below 
on two realizations: (i) = ff = f and (ii) = 0.54f and 

ff = 0.06f. (Since V is scaled in units of f, the trap confinement 
is effectively stronger in the second realization.) Next, we as¬ 
sume that the number of particles of each spin component No¬ 
can be tuned independently, as well as the interaction strength 
U between components can be controlled separately by means 
of Eeshbach resonances. In other words, one can always choose 
the optimal Hubbard parameters p and U from the conditions 
to approach and observe the largest AEM-ordered domains at 
a given initial temperature. According to the studies of two- 
component fermionic mixtures with hopping imbalance [2] we 
can set jUq = t//2, thus expect these domains in the trap cen¬ 
ter, and U = lOf* with t* = (ff + t[)l2 that leads to appear¬ 
ance of AEM correlations at high enough temperatures (up to 
T X 0.5f) in both cases. 

In Eig. 1 we show the dependence of the entropy per parti¬ 
cle S IN and, therefore, the initial temperature T of the trapped 



Figure 2: Real-space distributions of the on-site filling n, magnetization m = 
and the entropy per lattice site s in the harmonic trap with V = 0.015f 
obtained by LDA+DMFT in three specific points chosen from Fig. 1. 

Eermi gas before the lattice ramp on the temperature in the lat¬ 
tice T. We conclude that in mixtures with large hopping im¬ 
balance (the second realization) it is necessary to start with the 
lower initial temperature to approach the regime with AEM cor¬ 
relations. This may seem contradictory to expectations based 
on the entropy analysis of homogeneous systems [2]. How¬ 
ever, the effect becomes clear from the analysis of real-space 
distributions that are shown in Eig. 2. As one can notice, the 
effect originates from a smaller amount of entropy that can be 
distributed in mixtures with the hopping imbalance due to sig¬ 
nificant magnetization of metallic shells in the presence of the 
trapping potential. Note that, according to an additional anal¬ 
ysis performed and due to the real-space entropy distributions 
shown in Eig. 2, reducing the magnitude of the trapping po¬ 
tential in f/f* times and a corresponding increase of the total 
number of particles in the system does not help to suppress the 
main effect. 

Concerning a potential compensation of ferromagnetic 
shells, their presence is required to guarantee a proper param¬ 
eter regime in the bulk that is optimal for the easy-axis AEM 
state. Thus adding more particles of the less mobile (spin- 
down) component to the system will result not only in a com¬ 
petition between different types of AEM ordering in the bulk, 
but also in two ferromagnetic shells that are usually separated 
in real space (see Ref. [16] for details). The latter fact does not 
allow to significantly increase the entropy S jN necessary for 
observations of long-range AEM correlations. 

Therefore, inducing strong asymmetry in hopping ampli¬ 
tudes of two-component fermionic mixtures in a trap does not 
result in additional advantages from the point of view of ap¬ 
proaching the AEM-ordered regime with the highest possible 
initial temperature before the lattice ramp. However, a small 
(nonzero) asymmetry in hopping amplitudes may be still a very 
important ingredient for observations of the long-range AEM 
correlations as it is shown below. 

3.2. Direct observations of AFM correlations: importance of a 
specific direction for the symmetry breaking 

Let us consider the realization when the hopping amplitudes 
of both spin components can be considered as (a) approxi¬ 
mately or (b,c) strictly equal each other. In the first case, we 
consider that the difference in the hopping amplitudes is small 
not to affect thermodynamic properties of the mixture, so that 
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one can follow the obtained dependence for the balanced mix¬ 
ture shown in Fig. 1, but it is large enough to break the sym¬ 
metry in the Hamiltonian (1) in a controlled way. Accord¬ 
ing to Refs. [2, 16], at low temperatures the system prefers 
the easy-axis AFM-ordered configuration. In another limit of 
strictly equal hopping amplitudes we consider two realizations: 
(b) a small population imbalance is present, so that the system 
prefers the easy-plane AFM configuration, and (c) the mixture 
is balanced completely, so that the AFM state corresponds to 
spontaneous breaking of the SU(2) rotational symmetry in the 
pseudospin space. 

In the analysis we restrict ourselves to measurements based 
on QGMT, thus spatial correlations involving number of par¬ 
ticles of each spin component on all lattice sites that belong 
to a specific plane of the three-dimensional setup can be di¬ 
rectly observed. It means that the easy-axis correlations (in the 
Z direction of the three-dimensional pseudospin space) can be 
directly measured, but to obtain some signal from the easy- 
plane (XY-AFM) correlations one needs, e.g., an additional 
;7r/2-rotation in the pseudospin space to be performed before¬ 
hand (see also Ref. [6] for details). The latter in ultracold mix¬ 
tures of alkali-metal atoms is usually performed by applying 
additional rf-pulses. Therefore, from the experimental point of 
view, the above realizations of the symmetries correspond to (a) 
controlled both small nonzero hopping imbalance and the corre¬ 
sponding population difference (to produce weak ferromagnetic 
metallic shells in a trap), (b) zero (up to the experimental accu¬ 
racy) hopping and small nonzero population imbalance, and (c) 
balanced mixture in both respects. 

Now, to proceed with quantitative estimates in the framework 
of DMFT, let us specify the realizations (a) and (b). For the 
case (a) we take f-]- = 1.05f, T = 0.95f, and With 

the amplitude V - 0.015f this results in the polarization F = 
(Nf - Ni)/(Ni + Ni) 0.09 that remains almost constant with 
the temperature T jt (or, analogously, S fN) change in the region 
of interest. Hence, for the case (b) we take - ti - t and 
P ^ 0.09 that corresponds to A/j. - O.Olt. 

Next, by having access to the main quantity of interest, the 
on-site magnetization m, s as a function 

of the distance r to the trap center, we can reproduce the ex¬ 
pected single-site-resolution images for the occupancy by the 
single (e.g., spin-up) atomic component at specific values of 
S /N. To this end, we perform a Monte-Carlo (MC) sampling, 
where (m^ -H rijll) is used as a weighting factor for occupation 
of the lattice site i by the spin-up particle. In particular, by con¬ 
sidering two limiting cases at n - 1, (i) mf - +1/2, the site is 
always occupied (unoccupied) by the spin-up component, thus 
the checkerboard pattern is reproduced on the site /; (ii) mf = 0, 
the site i resembles the checkerboard structure or not with equal 
probabilities. 

Note that, while from the point of view of computational 
analysis and experimental observations the images for the case 
(a) can be obtained now in a straightforward way due to the 
easy-axis AFM order, in the case (b) the ordering takes place in 
the XY plane. Therefore, in the experiment one needs to per¬ 
form an additional 7t/2 rotation in the pseudospin space. At 
the same time, in the computational procedure, where the ro- 





Figure 3: Characteristic real-space single-shot images of the in-plane occupa¬ 
tion of lattice sites by one spin component at different entropies per particle 
and different experimental realizations (different types of the controlled sym¬ 
metry breaking) in a trap. The regions with n « 1, where AFM correlations 
can be present, are surrounded by circles. The images are taken close to the 
most probable values of F/Fq (see Fig. 4) from the corresponding sets of MC 
samples. 

tational t/(l) symmetry in the XY plane is broken to simplify 
the analysis (e.g., in the Xq direction), one needs to account 
for uncontrolled degrees of freedom (rotational symmetry) by 
multiplying the magnetization by a numerical factor that cor¬ 
responds to the averaging of the projection length of a random 
vector on the Xq axis, thus mf - Im^'' In. In the case (c) no ad¬ 
ditional rotations are necessary from the experimental point of 
view, but from the computational side the same arguments must 
be applied for the SU(2)-symmetric system, thus m^ - m^ jl, 
where Zq denotes the fixed direction in the numerical analysis. 

By considering all three cases (a)-(c) in Fig. 3 we directly 
see, that the effect produced by the controlled symmetry break¬ 
ing can be significant. In particular, we conclude that in the re¬ 
alization (a) at S IN - 0.6 the long-range correlations are much 
easier to observe directly by means of QGMT or Bragg scatter¬ 
ing. 

Finally, assuming that in the experiment one can perform 
an additional statistical analysis, in other words, one can av¬ 
erage over some finite set of real-space images (similar to those 
shown in Fig. 3) with the fixed system parameters (including 
temperature and the realization type), we can determine the av¬ 
erage size of the largest domain and analyze its dependence 
on the temperature and the type of the controlled symmetry 
breaking. For this purpose, we apply a procedure similar to the 
Hoshen-Kopelman algorithm in percolation theory [22] with 
accounting for the checkerboard structure of domains and sum¬ 
marize our results in Fig. 4. 
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Figure 4: Dependencies of the mean size of the largest in-plane AFM domain 
inside the region r < I2a (with the total number of sites Fq = 441) on the 
entropy per particle in the system. Each point corresponds to averaging over 
100 MC samples (real-space images). 

We conclude that this statistical analysis helps to extract 
more useful information about the AFM correlations in addi¬ 
tion to the single-shot images shown in Fig. 3. In all three re¬ 
alizations (a)-(c) the effects produced by magnetic correlations 
can be well identified and separated from the uncorrelated back¬ 
ground signal (with F/Fq ^ 0.13 in Fig. 4). However, even in 
this case we see a significant advantage of the realization (a) 
that has the largest signal-to-noise ratio in the whole region of 
entropies below S/N ^ 0.9^. 

4. Conclusions 

We studied theoretically perspectives of state-dependent op¬ 
tical lattices realized in Ref. [1] for the purpose of approaching 
and observing many-body quantum states with long-range mag¬ 
netic correlations in ultracold two-component fermionic mix¬ 
tures. We conclude that with account of the external harmonic 
potential there are no additional advantages of mixtures with 
large imbalance in hopping amplitudes for approaching quan¬ 
tum magnetism at the fixed value of the entropy per particle in 
the system. The strongest effect in this case originates from the 
structure of metallic shells that have nonzero magnetization and 
thus contain less entropy than in the balanced case. 

We analyzed a possibility to control the symmetry-breaking 
direction by using state-dependent optical lattices. It is shown 
that systems with small nonzero hopping imbalance have clear 
advantages among others for observations based on measure¬ 
ments of the particle number and, in particular, in the developed 
experimental techniques involving quantum gas microscope for 
ultracold fermions [7, 8, 9, 10, 11, 12]. We obtained charac¬ 
teristic single-shot images and performed statistical analysis of 


^Additional estimates based on DMFT analysis show that for a two- 
dimensional system with square lattice geometry the obtained dependencies 
remain qualitatively the same with a corresponding decrease in values of S /N 
(depending on the parameters of a particular experimental setup) approximately 
by 20%. 


their structure that allowed us to determine characteristic val¬ 
ues for the entropy per particle corresponding to appearance of 
sizable AFM-ordered domains in the system. 

The dependence on the smallness of asymmetries required 
for a controlled symmetry breaking can influence nonequilib¬ 
rium properties such as thermalization time for approaching a 
particular magnetically-ordered state. In this paper, by focusing 
only on equilibrium characteristics, we assumed that the intro¬ 
duced magnitudes of asymmetries in hopping amplitudes and 
number of particles are large enough for the system to equili¬ 
brate within the experimental timescale. Therefore, nonequilib¬ 
rium properties from that perspective remain highly important 
and could be studied in detail with further extensions [23] of 
the used theoretical approach. 
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